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ABSTRACT 

We present near-infrared spectroscopy of low-mass companions in a nearby triple system 
HD 130948 (Gliese 564, HR 5534). Adaptive optics on the Subaru Telescope allowed spectroscopy 
of the individual components of the (X'13 binary system. Based on a direct comparison with a 
series of template spectra, we determined the spectral types of HD 130948B and C to be L4 ± 
1. If we take the young age of the primary star into account (0.3-0.8 Gyr), HD 130948B and C 
most likely are a binary brown dwarf system. 

Subject headings: stars: low-mass, brown — stars: binaries: close — stars: individual (HD 130948) 



1. Introduction 

HD 130948 is a triple system discovered by Pot- 
ter et al. (2002) using the University of Hawaii 
adaptive optics system (Hokupa'a) in the course 
of an imaging survey of low-mass companions to 
nearby young stars using the Gemini North tele- 
scope. The HD 130948 (Gliese 564, HR 5534) 

1 Based on data collected at Subaru Telescope, which 
is operated by the National Astronomical Observatory of 
Japan. 

2 Visiting astronomer at the Institute for Astronomy, 
University of Hawaii 



system consists of a G2V primary star and a bi- 
nary companion at a separation of 2'.'63 from the 
primary. The binary companions are separated 
by 0'.'134, or 2.4 AU at a distance of 17.9 pc de- 
termined by Hipparcos. Potter et al. (2002) con- 
firmed the common motion of the members of the 
triple system with a baseline of 7 months to reject 
the chance projection of background sources. By 
using the AO system on the Subaru Telescope, we 
obtained individual spectra of each component to 
confirm the substellar nature of the faint compan- 
ions, HD 130948B and C. 
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2. Observations and Data Reduction 

The spectroscopic observation was made on UT 
2001 May 3 using the Infrared Camera and Spec- 
trograph (IRCS; Tokunaga et al. 1998; Kobayashi 
et al. 2000) with the 8.2 m Subaru Telescope in 
conjunction with its adaptive optics (AO) system 
(Takami et al. 1998; Gaessler et al. 2002). The 
Subaru AO system, which is equipped with a 36- 
element curvature sensor, is installed at the front 
end of the telescope Cassegrain port. A medium- 
resolution grism was used with a O'.'IO slit in the 22 
mas camera section of the IRCS to provide spectra 
with a resolving power of 800-1000 in the H and 
K bands. HD 130948A was used as the reference 
source for the AO system. The slit was put along 
the position angle of HD 130948B and C to cover 
them simultaneously. The spectra were recorded 
by nodding the telescope by 1" along the slit to 
subtract the sky emission. The total on-source in- 
tegration time was 1200 s in both the H and K 
bands. A nearby B8V star HD 164352 was ob- 
served as a spectroscopic standard at similar air- 
mass. The standard star is bright enough to give 
nearly the same AO correction as for HD 130948A. 
The spectroscopic fiat field was obtained at the 
end of the night with a halogen lamp. The seeing 
was 0'.' 4-075 at 2.2 fim throughout the observing 
period. 

Because of the close proximity and the large 
difference in luminosity between the primary and 
companions, the sky background in the spectrum 
is affected by the gradient in the wing of the point- 
spread function (PSF) of the bright primary. After 
a pair subtraction of the sky the residual back- 
ground slope was removed by interpolating the 
background on both sides of the companion. 

We obtained the one-dimensional spectra us- 
ing the IRAF 3 aperture extraction package after 
flat-fielding and bad-pixel correction. Special at- 
tention was paid to the adjustment of the aper- 
ture width. Figure 1 shows crosscuts of the ob- 
served spectrogram in the spatial direction with 
an overlay of the standard star profile normalized 
at the peak of each object. The FWHM of a stan- 
dard star is as sharp as 071 and 0708 in the H 

3 IRAF is distributed by the National Optical Astronomy 
Observatories, which arc operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 



and K bands, respectively. However, the inten- 
sity of the first diffraction ring is significant. The 
separation between the centroid and the peak of 
the first diffraction ring for an ideal diffraction- 
limited image of the 8.2 m circular opening should 
be 0707 and 0709 in the H and K bands, respec- 
tively. This is similar to the angular separation of 
0713 between HD 130948B and C. Consequently, 
the energy of HD 130948B in the first ring as well 
as the halo overlaps with HD 130948C. This may 
produce similar spectral signatures in both com- 
ponents that may be misleading. 

We calculated the contamination expected with 
the standard star PSF. The flux contamination of 
C by B within the 3 pixel (07067) aperture width 
is 25% and 22% in the H and K bands, respec- 
tively. The contamination of B by C is <20%. The 
extraction aperture width of the standard star is 
set to the same as that of the objects. To pro- 
vide the correct weight to the optimal extraction 
method (Home 1986), the statistical error in the 
2D spectrogram image is calculated by combining 
the readout noise and the shot noise from the sky 
emission and PSF spillover from the primary. The 
wavelength calibration was performed by fitting 
about 20 emission lines of the argon arc lamp of 
the Subaru Telescope calibration unit with a linear 
function. 

The intrinsic stellar lines in the standard star 
spectra, which are mostly atomic hydrogen lines, 
were difficult to remove because of their super- 
position with the telluric absorption lines. Only 
six of them at 1.588 ^m (BrlO), 1.641 fim (Br8), 
1.681 ^m (Br7), 1.737 /im (Br6), 2.166 fim (Br 7 ), 
and 2.374 /im (Pf20) were fit with a Lorentzian 
profile and subtracted before division. The small 
discrepancy in the airmass between the standard 
and the object was corrected by rescaling the stan- 
dard star spectra according to Beer's law. Rough 
flux calibration relative to the standard star was 
made in each band by assuming the intrinsic spec- 
tra of a B8V star is represented by a Planck func- 
tion of T cS = 10700 K (Tokunaga 2000). 

3. Result and Discussion 

3.1. Spectral Type of HD 130948B and C 

To determine the spectral types of HD 130948B 
and C, we compared the object spectra with tem- 
plate spectra of low-mass stars published by Reid 
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ct al. (2001), Geballe et al. (2002), Leggett et al. 
(2000), and Leggett et al. (2001). We found none 
of the template spectra match well with the object 
spectra. The continua of the object spectra were 
too steep to be a low-mass star, though the low- 
mass nature is apparent in the strong absorption 
of H2O at both ends of the H and K bands. We 
concluded that the continuum slope was affected 
by the AO correction. Not only a little offset from 
a narrow slit, but also slight difference in seeing 
can result in the variation of the continuum slope 
because of the sensitive wavelength dependency of 
the encircled energy on the AO performance (M. 
Goto et al., in preparation). The variation of the 
continuum slope from exposure to exposure was 
sometimes as much as 10% per 0.1 /im. Most of 
the time the difference was well represented by a 
simple linear function. 

Because of the uncertainty of the continuum 
slope, we had to correct the object spectra by mul- 
tiplying a linear function before comparison with 
templates. The strength of H2O absorption in- 
creases with spectral type, and the spectral change 
between any two spectral types is not a linear func- 
tion. In other words, we cannot convert the spec- 
trum of an LI star to an L8 star by just multi- 
plying a linear function. Therefore, in principle, 
we can determine the spectral type of an object 
uniquely by matching both the H 2 absorption 
with a template spectrum, even if we have a lin- 
ear continuum slope uncertainty. 

First, we assumed HD 130948B (or C) is a cer- 
tain spectral type, and estimated how much lin- 
ear continuum correction is required to best match 
the object spectrum with the template spectrum 
of the assumed spectral type. We calculated 
the cross-correlation between the object spectrum 
and the template spectrum as a function of the 
steepness of the applied linear function correction. 
The slope of the function is determined so that 
the cross-correlation between the corrected object 
spectrum and the template spectrum is highest. 
We repeated this procedure for each spectral type 
from LI to L8. 

Second, we investigated how well it matched 
each spectral type from LI to L8. The value of the 
cross correlation calculated above peaks at L3-L5 
in the H band and at L1-L5 in the K band. Con- 
sidering the simultaneous match in the H and K 
bands, we conclude both HD 130498B and C agree 



best with the template for L4 ± 1. This is consis- 
tent with the spectral type of L2 ± 2 determined 
independently by Potter et al. (2002). 

The final object spectra after correction are 
shown in Figure 2 with overlay of the best- 
matched template spectra, 2MASSW J0036+18 
(L4). The spectra were binned by 2 pixels 
along the dispersion direction. The spectra of 
HD 130948B and C turn out to be very similar 
to each other in detail and are almost twins, ex- 
cept that HD 130948B is 0.2-0.3 mag brighter 
than C. The deep absorption bandhead of the CO 
molecule at 2.3-2.4 fim is conspicuous. The pres- 
ence of absorption lines of Na I (2.209 /im) and 
K I (1.516 /zm) supports the idea that these are 
low mass stars. The absorption band of CH4 is 
not seen. This is in good agreement with the con- 
clusion that HD 130948B and C are not later than 
mid L-type stars. A series of the "unidentified" 
absorption features reported by Reid et al. (2001) 
at 1.58, 1.613, and 1.627 ^m can be all attributed 
to the FeH molecule (Wallace & Hinkle 2001). 
Ca I lines noted by Jones et al. (1994) and Tin- 
ney, Mould, & Reid (1993) arc marked in the 
figures. 

3.2. Mass of HD 130948B and C 

The range of spectral types derived above cor- 
responds to an effective temperature of T c g = 
1900 ± 75 (Leggett et al. 2001). On the basis of 
moderate X-ray activity, a fast rotation period (~ 
7.8 days), and photosphcric lithium abundance, 
the age of the primary should be 0.6 ± 0.2 Gyr 
(Gaidos 1998; Gaidos, Henry, & Henry 2000). 
On the other hand K. Fuhrmann (2001, private 
communication) finds that HD 130948 A may be 
a member of the Ursa Major Stream, which has 
estimated age of 0.3 to 0.65 Gyr (Giannuzzi 1979; 
Soderblom & Mayor 1993; Palous & Hauck 1986; 
Chen et al. 1997). 

Figure 3 is a reproduction of the evolutionary 
tracks of various low-mass stars of solar metallic- 
ity calculated by Baraffe et al. (1998). Regard- 
less of the uncertainty in the age determination 
of the HD 130948 system, the estimated masses 
for HD 130948B and C are under the sustainable 
hydrogen burning limit, 0.075 M for solar metal- 
licity, for both companions, and the most likely 
mass of the object is 0.040-0.065 M Q . 
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4. Summary 

We have presented moderate-resolution near- 
infrared spectra of the close binary companions 
recently discovered near HD 130948. The general 
spectral features of the companions, HD 130948B 
and C, are very similar. Based on a direct compar- 
ison with template spectra of low-mass stars, we 
have determined the spectral type of HD 130948B 
and C is L4 ± 1. This is consistent with the result 
of Potter et al. (2002), who found a spectral type 
of L2 ± 2 for the binary companions. If we take 
the young age of the primary star (0.3-0.8 Gyr) 
into account, then HD 130948B and C are binary 
brown dwarfs. 

We thank T. Geballe, S. Leggett, I. Ncill Rcid, 
and their collaborators for making the template 
spectra available in machine readable form for use 
in Figure 2. We also thank I. Baraffe and collab- 
orators for making the tabular data available at 
their ftp site for use in Figure 3. We thank all the 
staff and crew of the Subaru Telescope and NAOJ 
for their valuable assistance obtaining these data 
and their continuous support for IRCS construc- 
tion. M. Goto is supported by a Japan Society for 
the Promotion of Science fellowship. 
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Fig. 1. — Cross-section of the spectrogram across 
the slit length in the H and K bands with an 
overlay showing the spatial profile of the standard 
star. The standard star profile is registered to 
the position of each component and normalized at 
the peak of the object profile. The shaded areas 
represent the size of the aperture that the one- 
dimensional spectra are extracted from. 

Fig. 2.— H- and if -band spectra of HD 130948B 
and C normalized to the flux of HD 130948B at 
1.65 and 2.20 /im, respectively. The overlay is 
the best-matched template spectrum, 2MASSW 
J0036+18 (L4). 

Fig. 3. — The tracks of effective temperature for 
low-mass stars as a function of age (Baraffe et 
al. 1998). The effective temperature range for 
HD 130948B and C is indicated by the rectan- 
gular enclosure, along with the age determination 
uncertainty. 
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